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ABSTRACT 

Following on from our discovery of a significant population of M31 outer halo globular 
clusters (GCs), and updates to the Revised Bologna Catalogue of M31 GCs, we inves- 
tigate the GC system of M31 out to an unprecedented radius (« 120 kpc). We derive 
various ensemble properties, including the magnitude, colour and metallicity distribu- 
tions, as well as the GC number density profile. One of our most significant findings 
is evidence for a flattening in the radial GC number density profile in the outer halo. 
Intriguingly, this occurs at a galactocentric radius of 2 degrees (~ 30 kpc) which is 
the radius at which the underlying stellar halo surface density has also been shown to 
flatten. The GCs which lie beyond this radius are remarkably uniform in terms of their 
blue (V— I)o colours, consistent with them belonging to an ancient population with 
little to no metallicity gradient. Structural parameters are also derived for a sample of 
13 newly-discovered extended clusters (ECs) and we flnd the lowest luminosity ECs 
have magnitudes and sizes similar to Palomar-type GCs in the Milky Way halo. We 
argue that our findings provide strong support for a scenario in which a significant 
fraction of the outer halo GC population of M31 has been accreted. 

Key words: galaxies: star clusters - galaxies: interactions - galaxies: formation - 
galaxies: evolution - galaxies: individual (M31) - galaxies: haloes 



1 INTRODUCTION 

The properties of globular cluster (GC) systems provide 
valuable prob es of the formation and evolution of their host 
galaxies (e.g. IWest eralll2004l : ferodie fc Stradeij|2006l ). It 
is commonly believed that GCs form in major star-forming 
episodes that accompany g alaxy formation, as wel l as in sub- 
sequent merger events (e.g. IZepf fc AshmanlllQQ^ ). Further- 
more, the native GC population of a galaxy is expected to be 
augmented through mergers with and accretions of smaller 
systems, each of which will bring its own retinue of GCs into 
the final galaxy. As a result, the GC population of a galaxy 
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will reflect both the amount of mass formed in situ, as well 
as that which has been accreted. 

As GCs are (mostly) luminous and compact, they can 
be readily observed in galaxies up to a few hundred Mpc dis- 
tant. Nevertheless, the GC systems of Local Group galaxies 
remain of central importance as they allow the most de- 
tailed study of the properties of GC populations and how 
they correlate with the formation history of their host galax- 
ies. This is possible as their proximity permits studies of 
resolved field and cluster stellar popul ations th r oug h their 
colour-magnitude diagrams (e.g. iMackev et al.l l2006h and 

spectroscopy (e.g lBarmbv et aLl2000t ). ' 

The study of the Galactic GC system bv lSearle fc ZinnI 

(|l978t ) proved crucial for understanding the history of our 
own Milky Way (MW). Their analysis of GC metallicities 
within the halo led them to conclude that many of these 
objects must have formed within protogalactic fragments 
that fell into the Galaxy after the collapse of the central 
regions had been completed. This scenario was in stark con- 
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trast to the monolithic slow collapse picture earlier pro- 
posed bv lEggen. Lvnden-Bell. fc Sandagj (|l962l ). One of the 
main distinguishing characteristics of the two models was 
the halo metallicity gradient; the latter scenario predicted 
a radial gradient with clusters nearer to the centre being 
more metal-rich, as they were formed somewhat later in the 
collapse, while the former predicted a spread in metallic- 
ities at all radii but no significant radial variation. Mod- 
ern theories of galaxy assembly are based on hierarchical 
structure form ation with galaxies fo rming inside dark mat- 
ter halos (e.g. IWhite fc Frenklll99ll ). Ahhough this model 
bears some similarities to the scenario proposed by Searle & 
Zinn, there are also several differences. In particular, it has 
been shown that build-up via accretion can sometimes lead 
to halo metallicity gradients since massive satellites, which 
are normally more metal-rich, sink further into t he potential 
well of the host than do low mass objects (e.g. I Font et aU 
l2008l : lDe Lucia fc Helml2008 ). Thus, while the lack of a halo 
metallicity gradient argues against slow pressure-supported 
collapse, the existence of one could be consistent with both 
that scenario as well as accretion. 

The GC system of M31 has also been studied intensively 
in order to searc h for clues about how that system formed 
and evolved (e.g. Crampton et al.l 19851: Elson fc Walterbo3 



19881: iHuchra. Brodie. fc Kent! Il99ll : iBarmbv et al.l l20o'oF 
Perrett et al.ll2002l : iFan et al.ll2008f rSince M31 is similar to 
the MW in many respects, it may be expected to have expe- 
rienced a similar assembly history. Early work suggested a 
mild gradient in the metallicity of the M31 GC system (e.g. 
Sharovlll988l: iHuchra. Brodie. fc Kentlll99ll: [Sarmbv et al.l 



2000l . IPerrett et al.l (|2002l ') and lFan et al.l (|2008''l found that 
this result is primarily driven b y a metallicity grad ient in 
the metal-poor GCs alone, while IPerrett et aU l|2002r ) noted 
that the slope of the gradient may flatten beyond a projected 
radius of ~ 14 kpc. 

In addition to the metallicity gradient, another key 
property of a GC system is the radial profile - the areal num- 
ber density of GCs as a function of distance from the cen- 
tre of the host galaxy. Radial GC profiles have traditionally 
been represented with either a power-law (i.e. R~") or a de 
Vaucouleurs law (i.e. R^^^ law), with a flattening of the pro- 
file at small radii, and a gradu al steepening in the outer re- 
gions l|Brodie fc Stradeill2006l ). The most re cently published 
GC su rface density profile in M31 is that of iBattistini et ahl 
l|l993ll . who experimented with a variety of fitting formulae. 
Using a sample which extended to a radius of ~ 30 kpc, they 
noted that their data were consis tent with a st eepening in 
the outer regions, earlier found bv lRacinell|l99ll ). if one uses 
the general power - or R^^* laws (although this is not the 
case if they employ an R^^^'^ law). They argued that this 
steepening was unlikely to be due to incompleteness in their 
sample at large radius, although the possibility could not be 
completely excluded. 

A major limitation of all previous studies of the ensem- 
ble properties of the M31 GC system has been the restricted 
radial range of the samples. These studies have generally em- 
ployed samples extending to no more than ~ 25 kpc, while 
analogous studies of the MW GC system have extended to 
beyond 100 kpc. In the last decade, vast amounts of new 
imaging have been obtained of the outer regions of M31 re- 
vealing that both fi eld stars and GCs extend to radii of well 
over 100 kpc (e.g. iFerguson etHI I2OO2I : llbata et al.l l2007l : 
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Figure 1. The location of the new globular clusters (red trian- 
gles) and extended clusters (blue stars) in relation to confirmed 
RBC GCs that lie in the survey area analysed here (black cir- 
cles). The limited extent of the known GC population prior to 
this study can be clearly seen. The dashed line outlines the INT 
survey area covered, and the dashed-dotted line outlines that part 
of the Megacam survey employed here. The inner ellipse has a 
semimajor axis of 2° (27 kpc) representing a disk with an in- 
clination of 77.5°; the optical disk of M31 lies well within this 
boundary. The outer ellipse denotes a flattened ellipsoid of semi- 
major axis length 4° (55 kpc). The dashed circle lies at a radius 
of 30 kpc, and shows the region at which the break in the GC 
surface densi ty profile occurs (see text for details). The outermost 
GC fou nd bv lMartin et al.l 1I2OO6I') and discussed bv lMackev et all 
||201o3) is also shown (M06 GC) and indicates the extent of the 
currently-known M31 halo GC population. 



Huxor et all I2OO8I : iMcConnachie et all l2009l : iMackev etahl 
2010al 'l. It is therefore timely to revisit the properties of the 
M31 GC system in light of these new datasets. In Section 
2, we discuss the sample used for our study which consists 
of the revised Bologn a catalogue and new GCs presented 
in iHuxor et al.l l|2008l) (hereafter Paper I). In Section 3, we 
discuss radially-dependent properties of the M31 GC system 
using a baseline that extends over 100 kpc. In Section 4, we 
discuss properties of extended star clusters in M31's outer 
halo. Section 5 discusses our findings in the context of the 
assembly history of M31. 



2 THE SAMPLE 

The sample of GCs used in this study is taken from ver- 
sion V3.33 of the Revised Bolo gna Catalogue (RBC), which 
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was released in March 20ojf]. A fu l l des cription of the 
R BC can be found i n iGalleti et al.l (12004) w ith updates 
in'Galle ti et all ^0&) and lGalleti et all l|2007l ). It contains 
1983 entries, including 509 confirmed GCs, 13 confirmed ex- 
tended clusters (ECs) and 1049 "candidate clusters". This 
ve rsion of the RBC i ncludes 103 confirmed objects reported 
bv lKim et all (|2007l ') although subsequent work has shown 
that many of these objects have been imp roperly classified 
l|Caldwell et al.ll2009l : [Peacock et al.ll2010l V This has only a 
minor impact on the results presented here as the Kim et 
al. objects all lie within a projected radius of 18 kpc (90% 
within 10 kpc) from the centre of M31, while our main fo- 
cus in this paper is on the halo regions beyond 20 kpc. The 
RBC V3.5 also includes the 40 new outer halo clusters found 
by our group and reported in Paper I. In brief, these lat- 
ter objects were discovered in more than 80 square degrees 
of imaging survey data taken with the Wide-Field Cam- 
era on the Isaac Newton Telescope and Megacam on the 
Can ada France Hawa ii Telescope (see iFerguson et al]|2002l . 
and llbata et akllioOTl . for details of the surveys). Thanks to 
the quality of the seeing conditions and the proximity of 
M31 to us, the clusters could all be identified as marginally 
or fully resolved stellar concentrations. The new clusters we 
discovered increased the total number of confirmed GCs in 
M31 known at the time of publication by a modest amount 
(~ 10%) but the number of confirmed GCs known beyond 
1° (f» 14 kpc) by more than 75%. 

For the main sample used in our analysis, we select clus- 
ters in the RBC V3.5 that have a Global classification flag 
(f) equal to 1 or 8 corresponding to 'confirmed' GCs and ECs 
respectively. Although the RBC uses this flag to distinguish 
between ECs and GCs, studies to date sug gest M31 ECs 
have the same ancient stellar content as GCs ( Mackev et al.l 
I2OO6I ). hence we consider them to gether in the p r esent work. 



The very remote GC found by Martin et al. (|2006|) and 



recently studied in detail by iMackev et al. ~(|2010al ) 



eluded in the sample for most of the analysis but not for 
the GC surface density profile as it lies in an area of the 
CFHT/Megacam survey for which full results of our GC 
search have yet to be published. Possible young GCs are also 
excised from our sample by removing those clusters for which 
the young cluster flag (yy) is greater than zero. While M31 
appears to have a g enuine population of young to interme - 
diate age GCs (e.g. iFusi Pecci et allliooa iFan et al.ll200d ). 
such a population is absent in the MW. Since our primary 
interest in this paper is to examine ancient star clusters as 
probes of galaxy formation, a sample with these young ob- 
jects removed makes for a more appropriate comparison of 
the GC systems of both galaxies. The final sample we em- 
ploy for the bulk of our analysis contains 431 objects while 
the sample used for constructing the surface density profile 
analys is contains 430 objects - removing the iMartin et al.l 
l|2006t ) GC as noted above. 

The RBC is a compilation of a variety of different data 
sources and in situations where our own photometry is avail- 



^ While this paper was in preparation, the RBC was updated in 
December 2009, to V4.0. However the changes do not affect the 
results in the present work, as the updates concern a few GCs in 
the inner region of M31 and/or add new young clusters. We do 
not include such young clusters from our sample (see main text). 



able from iHuxor et al] (|2008l ). this is used in preference to 
that provided by the RBC for homogeneity. The projected 
galactocentric radii of the GCs are re-derived from the equa- 
torial coordinates listed in the RBC using M31 central co- 
ordinates of RA= 00''42'"44.3'' and Dec= -h41°16'09" for 
the centre of M31, taken from NEE0. Figure [T] shows the 
spatial distribution of the GCs around M31 in our sample 
and illustrates the much greater radial extent of our sample 
compared to those used for previous M31 GC studies. 



3 RADIALLY-DEPENDENT PROPERTIES OF 
THE M31 GC SYSTEM 

We now proceed to use this sample to derive the basic 
properties of the M31 GC system out to much larger ra- 
dius than has been previously done. Wherever possible, we 
compare the quantities derived with those of the MW GC 
system for which we use the information listed in t he Mc- 
Master Catalogue of Milky Way Globular Clusters Q (|Harrij 
This catalogue lists data for 141 GCs, although not 
all of them have a full set of derived parameters. The cat- 
alogue also does not contain a few very recently-discovered 
MW GCs in the G alactic plane (e.g. lKobulnickv et al.ll2005l : 
iKurtev et al.|[200^ ) but since these do not lie at large galac- 
tocentric distances, they do not affect our overall results. We 
do, however, include the t wo newly-discovered G Cs found in 
the outer halo of the MW (|Koposov et al.ll2007l ). As there is 
limited information available for these objects at the present 
time (only size and V-band magnitudes), they can only be 
used for a subset of the following comparisons. The values 
of E(B-V) given in the Harris catalogue are used to derive 
the corrected magnitudes and colours for each MW cluster. 



3.1 GC Luminosities 

The M31 GC luminosity function (GCLF) is shown in Fig- 
ure[2l The M31 GC magnitudes have be en calculated assum- 
ing a distance modulus of 24.47 mags (|McConnachie et al.l 
l2005h . The magnitud es are corrected for ex t inctio n using 
the E(B-V) values in iFan. de Griis. fc Zhoul (|2010l ) where 
available. These values are derived from SED fits to multi- 
band photometry and give the total line-of-sight extinction 
including that internal to M31. For those GCs in our sa mple 
that were not studied bv lFan. de Griis. fc Zhoul (|2010l ). we 
adopt their median value of E(B-V) = 0.12, except for the 
very outer halo GCs (Rproj >40 kpc), where we use only 
Galact ic foreground extinction calcu l ated b y interpolation 
of the ISchlegel. Finkbeiner. fc David (Il998l) mapiQ. Previ - 
ous work on the M31 outer halo GCs (|Mackev et all [20071 ) 
gives reddening values consistent with Galactic extinction 
being the significant component. The MW GCs are overlaid 
in red. 

Although the number of currently-confirmed M31 GCs 
is already ~ 3 times larger than the number known in the 
MW, the shapes of their luminosity functions are roughly 
similar, at least down to an absolute magnitude of Mvo ~ 



^ http:/ /nedwww. ipac.caltech.edu/ 

^ http:/ /physwww. physics. mcmaster.ca/f^harris/mwgc.dat 
^ http:/ /astro. berkeley.edu/~marc/dust/data/data. html 
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Figure 2. Extinction-corrected absolute V-band magnitudes 
of M31 GCs (black line) in our catalogue, using val- 
ues of Ay derived from IPan. de Griis. &: Zhoul | |201C| ) or 
ISchlegel. Finkbeiner. &: Davi depending on availability 

(see text for details). The red dashed histo gram shows the distri- 
bution for the MW GCs, and includes the lKoposov et alj ||2007| ) 
clusters. Median values for each sample are shown with vertical 
lines (-7.9 for M31 and -7.3 for the MW). 



—5, where completeness issues begin to complicate the M31 
sample. Figure [5] reveals an offset of ~ 0.6 mags between 
the peaks of the GCLFs (the medians being -7.9 and -7.3 
for M31 and the MW respectively) however the magnitude 
of this offset is very dependent on the assumed extinction 
for th e M31 population. While the ,Fan . de Griis, fc Zh ou 
reddening values are the best available as they are 
derived in a homogeneous manner for a large number of clus- 
ters, there are hints that they may overestimate the extinc- 
tion at high E(B-V) when compared to other methods (see 
their Figure 7) . Indeed, if only Galactic foreground extinc- 
tion is assumed, there is no discernible offset in the GCLF 
peaks of M31 and the MW. The M31 GCLF also exhibits a 
secondary peak at Mvo ~ —5.5 that is not seen in the MW 
population. This peak is present when inner (^ 20 kpc) and 
outer (^ 20 kpc) halo samples are considered separately, 
sug gesting it c o uld b e genuine. On the other hand, many of 
the lKim et al] l|2007l ) clusters have magnitudes around this 
value, a large fraction of which have questionable cla ssifi- 
cations (e.g. see the discussion in iPeacock et al.1 (|2010l )). A 
more detailed comparison of the MW and M31 GCLFs will 
be carried out at a later date, when uniform photometry 
and reddenings are available for the entire sample and the 
completeness and contamination of the M31 catalogue has 
been more rigorously quantified. 

Figure [3] shows the distribution of absolute magni- 
tudes against projected galactocentric distance for those 
GCs for which V-band data is available. This includes all 
but four of the Harris MW catalogue and all but eight of 
the M31 catalogue. The galactocentric distances for the MW 
GCs, Rgc, are converted to equivalent projected radii using 
roj — 7r/4 X Rgc- This numerical factor accounts for the 
ratio of the projected distance to the true distance for an 



Figure 3. The radial variation of the extinction-corrected V-band 
absolute magnitudes for M31 (open diamonds) and MW (inverted 
triangles and open squares) GCs, using the same data as Fig. (2] 
In the case of the MW GCs the actual distance (Rgc) is converted 
to an "average projected distance" through Rproj = Rgc x {tt/A). 
The unusual MW cluster NGC 2419 is labeled. 

isotropically-distributed GC population viewed from a dis- 
tant random external vantage point. 

Figure [3] reveals that M31 has a population of lumi- 
nous GCs at large galactocentric distances, a finding pre- 
viously co mmented on from studies with smaller datasets 
iMackev et al. 20oi [Calleti et al. 2001). While these GCs 
are no more luminous than GCs at smaller radii in M31, they 
are significantly more luminous than the outer halo GCs in 
the MW. In the MW, only NGC 2419 has a galactocentric 
radius larger than ~ 35 kpc, and a luminosity greater than 
Mvo <j -6.5 whereas M31 has 21 GCs in this same region 
of parameter space (even although the sample studied here 
covers only a quarter of the sky area at these distances). Al- 
most all the outer halo GCs in the MW are "Palomar-type" 
clusters with low luminosities and relatively large half-light 
radii. As both samples are likely to be highly complete in 
this magnitude range, the difference between the numbers 
of luminous M31 and MW GCs, a factor of ~ 20, cannot 
be solely due to the difference in the overall size of the GC 
populations (which accounts for only a a factor of ~3). It is 
unclear whether M31 also possesses an excess of faint outer 
halo GCs compared to the MW since completeness issues 
currently complicate our analysis below Mvo ~ -5. 

3.2 GC Colours and Metallicities 

Fig |4] shows histograms of the (V-I)o colours for the MW 
and M31 GC samples. Only 96 out of the 141 MW GCs had 
(V-I) values listed in the Harris database. As an additional 
20 MW GCs had (B-V) data, we fit the (V-I) colours as a 
function of (B-V) for those GCs for which both colours were 
available. This resulted in the relation (V-I) = 1.23 x (B- 
V) + 0.09 (RMS = 0.05), which was used to estimate (V-I) 
values for those remaining 20 GCs. Figure H] shows that the 
M31 GC system has the same median colour as that of the 
MW, with the reddening-corrected median (V-I)o colours 
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Figure 4. The distribution of (V-I)o for GCs in M31 (black 
solid histogram) and the MW (red dashed histogram), where this 
colour measure is available. Where no I band data are avail- 
able for MW GCs, (V-I) colours are derived from (B-V) (see 
text for details ) . Med ian values are shown by vertical lines. The 
IKoposov et al. I 1 I2OO7I ) GCs are not included here as no colour 
information is available for them. 



being 0.89 and 0.93 mags for M31 and the MW respectively. 
This is somewhat at odds with earlier findings that the M31 
GC system is, on average, redder than that of the MW 
(e.g. iHuchra. Brodie, &: Kentiri99ll ) and depends largely on 
our a dopted reddening values from iFan. de Griis. fc Zhoul 
(|2010l ). 

When viewed as a function of galactocentric radius (Fig- 
ure [5]), the mean GC colours in the outer halo are uni- 
formly blue and similar in both systems. We find ((V-I)o = 
0.87±0.04 in M31 for R> 30 kpc and (V-I)o = 0.92±0.03 for 
R> 15 kpc in the MW, where errors are the standard error 
on the mean. The colours are almost identical for 45 kpc, 
the difference for R> 30 kpc solely due to the inclusion of the 
very blue extended cluster HECl. As small-scale extinction 
variations are not likely to be a factor in these remote parts, 
this suggests that intrinsic properties of the outer halo GCs 
in both systems are very similar. Interestingly, the mean 
colour of the oute r halo GCs in M33 ((V-I)o = 0.88 ± 0.05 
iHuxor et al]|2009l is also very close to the values found in 
M31 and the MW. 

The dispersion in colour of the outer halo GCs in Fig- 
ure [5] is relatively small beyond 30 kpc and there is no ev- 
idence of any radial gradient. Since the integrated colour 
reflects both the age and the metallicity of a GC, further 
information is required in order to properly interpret the 
uniformity in the outer halo. For several of the M31 outer 
halo GCs, we have previously derived metallicities from fit- 
ting model isochrones and MW globular cl uster fiducials 
to high-quality co lour-magnitude diagrams (|Mackev et al.l 
I2006l . l2007l . l2010al 'l. These fits are based exclusively on the 
red giant and horizontal branch stars. Figure [6] shows the 
metallicities derived in this manner for 11 compact and 4 
extended GCs in the M31 halo. A linear least-squares fit 
to all the data yields a slope of -0.006 ± 0.002 dex kpc~^, 
indicating a marginal negative gradient. However there is 



Figure 5. The distribution of (V-I)o for M31(open diamonds) 
and MW (inverted triangles) with galactocentric radius. The data 
used here is the same as for Fig|4] The mean colours for the Kproj 
>30 kpc M31 and Kproj > 15 kpc MW GCs are shown as black 
and red dashed lines respectively. 

one clear outlier in the plot - H14 - which was previously 
fiagged by iMackev et al.l |20o3) as a possible intermediate- 
age cluster. If true, this would affect the metallicity of the 
cluster which has been derived assuming a 10 Gyr age. If 
this object is removed, we obtain the dashed line fit with 
a mean metallicity of -1.94 ± 0.22 dex and no discernible 
gradient (-0.003 ± 0.002 dex kpc^^). If we focus only on 
those GCs which lie beyond 30 kpc, the mean metallicity is 
effectively the same, given the errors, at -1.95 dex. It thus 
seems fair to assume that the uniformity in the broadband 
colours of the outer halo GCs seen in Fig.[5]reflects a narrow 
spread in both age and metallicity, with this population be- 
ing predominantly old and metal-poor. For comparison, the 
average [Fe/H] for the MW sam ple with Rproj > 15 kpc 
is -1.70±0.23 dex (jHarrisI [l99i ) . Our finding of no signifi- 
cant metallicity gradient in t he M31 halo GC popula tion is 
consistent with the study of I Alves- Brito et all (1200^), who 
found the same result from spectroscopic measurements of 
a small sample of GCs, although we note that there is a 
discrepancy between the absolute metallicities obtained by 
these authors and those inferred from our CMD fitting). 



3.3 The Radial Surface Density Profile of the 
M31 GC System 

The radial number density profile of the M31 GC system 
can be calculated in a straightforward fashion once a correc- 
tion is made for the non-uniform azimuthal coverage of the 
outer halo surveys from which our GC catalogue has been 
constructed. Indeed, as can be seen in Figure [1] the surveys 
used to identify outer halo GCs extend significantly further 
in the south eastern direction than anywhere else. Using 
knowledge of the field centres and sizes, we determined the 
proportion of the sky imaged within a specific annulus and 
hence the correction factor needed to estimate the number 
of clusters within that radial range (on the assumption of an 
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Figure 6. The radial variation of metallicities for t hose GCs with 
accurate values in t he literature from CMD fitting iMackev et alj 
||2006| . I2OO7I. Hoioi). The solid line is a fit to all the data, and 
the dashed line is a fit that exc ludes the o u tlier, the possibly 
intermediate-age cluster H14 ( .Mackev et alJ |20o3), where is it 
called GC7). The errors on indiv idual data points are ±0.15 dex 
jMackev et al.ll2006l. |2007|. l2010at) . 



isotropic distribution). Our number densities are calculated 
within circular annuli, which should be appropriate if the 
true halo shape is roughly spherical. 

Figure [7] shows the projected GC number density as a 
function of radius for M31, with the correction for incom- 
plete azimuthal coverage included. The data are presented 
in equally-spaced logarithmic bins. The ilcr errors are de- 
termined assuming the number of GCs at a given radius is 
a Poisson process. Then in each bin, the upper (lower) end 
of the error bar is set to the value of underlying GC density 
for which the probability of observing a number of clusters 
as lo w (high) as we actually did, or lower (higher), is 15.87% 
iMuld cr 1981)- 

The profile shows a broken power-law behaviour. Inside 
a projected radius of ~ 5 kpc, the radial number density 
profile of GCs is rather fiat, as h as been previously com- 
mented on in earlier studies (e.g. de Vaucouleurs fc Butal 
19781: [Harris fc RacindlToTgl : IWirth. Smarr. fc Brunclll985l : 
Battistini et al.l 19931 ). While this behaviour could be gen- 
uine, it could also represent incompleteness in the samples 
at small radii. Indeed, detecting GCs against the high sur- 
face brightness bulge which dominates in this region remains 
a challenge even for modern day surveys. Beyond this inner 
region, the profile is much steeper out to a projected radius 
of ^ 30 kpc after which it fiattens again. This is the first 
time that this behaviour has been seen in M31 and is due to 
the fact that the new profile extends to radii of ~ 100 kpc, 
m ore than three tim es further out than the earlier profiles 
of lRacind l|l99ll ) and lBattistini et all l|l993h . 

A broken power-law fit is overlaid in Figure[7]with best- 
fit indices of -0.80 ± 0.10 inside of ~ 5 kpc, -2.88 ± 0.13 
in the range 5 <^ Rproj <^ 30 kpc and —0.87 ± 0.52 outside 
of ~ 30 kpc. Although this power-law representation is typ- 
ically employed in extragalactic GC studies, it is clearly a 
rather poor description of the behaviour in M31, especially 



Figure 7. A log-log plot of the radial number density profile of 
GCs in M31, with Poisson errors (see text). A broken power-law 
has been fit, with slopes and errors for each component given in 
the legend. The inner region extends from the centre of M31 to 
a projected distance (Rproj) oi ^ b kpc (dotted black line); the 
intermediate region from Rproj ~ 5 to ?s 30 kpc (red dot-dashed 
line); and the outer region from Rproj ~ 30 kpc (solid black line). 
A change in the slope of the profile can be seen at 5 and 30 
kpc. The outer three bins contain 13, 10 and 3 GCs in the original 
data which increase to estimated values of ~ 15, 25 and 27 GCs 
when the incomplete spatial coverage for each annulus is taken 
into account. 



at intermediate radii, being too stee p at small radii and too 
shallow at large radii in this range. iBattistini et al.l (Il993l ) 
also noted that power-law fits were a poor representation 
of the radial number density profile of their sample of M31 
GCs and investigated a number of other empirical fits. They 
found most success with an R^''™ law, with m ~ 1.6. In Fig. 
[S] we show a profile of this form fit to our current sam- 
ple. Even though there have been many updates to the M31 
GC catalogue since Battistini's work, this form still provides 
an excellent description of the radial number density profile 
over a significant radial range (5 < Rproj < 30 kpc). Indeed, 
this relation even provides a good fit for the outer halo GC 
profile, albeit with a different slope. 

It is worth commenting that while the outer halo region 
is particularly rich in the more extended clusters (described 
below), these objects are not the primary cause of the flat- 
tening in the profile. Of the 36 GCs beyond the break in 
the GC surface density profile, only nine are ECs. This can 
be compared to 2 ECs from the 36 GCs in an annulus im- 
mediately interior to the break. If the radial surface density 
profile is constructed using only compact GCs, the the break 
is still clearly present, but the slope at large radii is some- 
what steeper. 



4 THE EXTENDED CLUSTER POPULATION 

One of the main results thus far from our GC search has been 
the discovery of a populati o n of E Cs within the halo of M31 
(Paper I and iHuxor et al.] (|2005l )). These objects typically 
have half-light radii of 20 pc and are significantly more 
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Figure 8. The radial profile of surface density of GCs in M31 
against R^^^'^, using the same bins as for Figure [7] Linear fits 
have again been made to the inner, intermediate and outer regions 
as in Fig. [Jl The fiattening in the last three bins remains very 
clear. 



extended than the no rmal GC populatio n w hich have half- 
light r adii of a few pc. lHuxor et al.l (|2005l ) and lMackev et al] 
l|2006t ) presented an analysis of the four brightest ECs dis- 
covered in the early stages of our survey. Here we present 
structural parameters for all thirteen ECs in the sample 
studied here. 

To investigate the structures of these clusters, we con- 
sider empirical King profiles. 



E(r-) = Eo 



l + {r/R,fy^^ {l + {Rt/R, 



where Rc and Rt are the core and tidal radii respectively. 
These profiles were fit to V band (in the case of INT) or 
g- band (for Megacam) photome try w h ich w ere the shortest 
wavelength data available to us. iKind (|l985l ) notes that the 
best cluster profiles are obtained by using the bluest images 
since statistical fluctuations due to individual red giants can 
be problematic at longer wavelengths. The fits were made 
to the cumulative fiux measured within a series of apertures 
of increasing radii made with the IRAF/apphot photometry 
packag^. Hence the fit, which minimised x^, was made to 
the integral of the King profile, which gives the total flux 
within a radius r: 



C{r) 



Eo27r 



R: 



2Ii: 



^ln(a)-^(a)i/^ + — + ^ 



/5 



2/32 



2RJ 

13 



where 

a = l + {r/R^f 
and 



P^{l + {Rt/R.fY^'' 

There are some uncertainties associated with the val- 
ues derived from our King profile fits. In the case of HEC7, 
contamination from foreground stars was removed first by 
using the "patch" option in to replace the affected 

region with an average of the background sky. Another clus- 
ter, HEC4, lies at the edge of a CCD chip but the fractional 
area missing is small allowing us to still obtain photome- 
try to a radius of 13 arcsec. Cluster HEC6 was so heavily 
affected by background galaxies that no proflle flt could be 
attempted. In this case, a value for Rh was obtained by flnd- 
ing the aperture containing half the luminosity as given by 
our measurement of the total magiutude within a 12 arcsec 
aperture, and thus has a signiflcant degree of uncertainty. 

For four of the ECs studied, we can compare the ground- 
based proflle flts made here with those derived from higher 
spatial resolution HST/ACS data and reported in Tanvir et 
al. (2011). The HST data allow for better removal of back- 
ground galaxy contamination and give proflle flts that are, 
for three of the four ECs, somewhat smaller than those found 
with the ground-based data. It is likely that this is mostly 
due to the limitations inherent in constructing proflles from 
ground-based data. This underscores the need to treat the 
quantities derived for faint ECs from ground-based flts with 
some caution. 

It is of interest to compare the size distribution of ECs 
with that of the more compact GCs in M31. Although we 
would have liked to obtain proflle flts for our entire sam- 
ple of outer halo GCs, this proved to be impossible for the 
compact objects. The average seeing of our ground-based 
data was 1.2 arcsec, which corresponds to a physical size 
of ~4-5 pc at the distance of M31. Thus, we do not re- 
solve classical GCs, which typically have core radii of 2-3 
pc, to a sufflcient extent to accurately flt their radial pro- 
flles. Instead, we use as a comparison sample the work by 
iBarmbv et al.l (2007) and Tanvir et al. (2011) which presents 
King proflle flts for compact M31 GCs from high resolu- 
tion imagery with HST. We also include the high resolution 
gr ound-based imaging st udy of the outermost M31 halo GC 
bylMackev et ajj (|2010al '). Figure [5] shows the half-light radii 
(Rh) distribution of M31 ECs (black (green) open diamonds 
show ground-based (HST) measurements) and GCs (crosses) 
against projected galactocentric radius compared to the MW 
GCs (flUed diamond s). We only show those members of the 
iBarmbv et al. I I2OO7I) sample that are also found in our mas- 
te r M31 GC cat a logue described previously. Most of those 
m IBarmbv et al] l|2007f) that are not in our list are putative 
"young" clusters, as indicated by the corresponding flag in 
the RBC database. A further handful are not yet included in 
the RBC. Figure [9] includes the newly found compact M31 
GCs from Paper I (plus signs) which are plotted with an up- 
per limit on the half-light radius of 4.5 pc. (black crosses), 
derived from the average seeing. The newly discovered GCs 
from Paper I with structural parameters derived from our 
HST data (Tanvir et al. 2011) are shown as green plus sym- 
bols). It can be seen that the M31 GCs exhibit a very differ- 
ent size distribution compared to those of the MW, in that 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



^ GAIA, the " Graphical Astronomy and Image Analysis" tool 
is now available as part of the Starlink Software Collection, via 
|http: //star link, ja ch . hawaii .edu/ 1 
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Table 1. Properties of the Extended Clusters. Columns 2-5 and 8-9 derived by fits to a King model. Columns 6 and 7 show values 
of Rc and Rh derived from recent HST/ACS imaging of HEC4,5,7 and 12, indicating errors in the ground-based data. The values for 
HEC6 are not derived from model-fitting, sec text for details. The cross identifications for the ECs in lHuxor et al. I ll2005l l and Tanvir et 
al. (submitted) are also given (H05 ID and TIO ID, respectively) where appropriate. 
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flhHST 
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Figure 9. The distribution of for M31 and MW GCs with pro- 
jected galactocentric distance (kpc) showing the MW GCs (filled 
dia monds), new ex tended clusters from Paper I (open diamonds) , 
the lBarmbv et al ] 12007) clusters that are in our catalogue (trian- 
gles). The new compact GCs from Paper I (plus signs) are shown 
with an value of 4.5 pc (except for ten clusters, where values 
are available from Tanvir et al. (submitted)). This is intended to 
show a maximum value, as they are unresolved in our data. Sec 
text for details. 



Although consistent with a globular star cluster, this result 
has large uncertainties and cannot be used to definitively 
exclude the presence of a modest amount of dark matter. 
The first three extended clusters reported in iHuxor et al.l 
l2005l were extreme in their large magnitudes and half-light 
radii. While these objects appeared to be somewhat iso- 
lated in the Mv-Rh plot, the additional clusters reported 
here indicate that ECs actually form a continuous, nearly 
vertical, sequence which overlaps at the faint end with the 
smaller, less-luminous Palomar-type GCs found in the MW. 
Moreover, the new HST data suggests that the most ex- 
treme ECs are not as large as previously thought. How- 
ever, they are still significantly more extended than typi- 
cal GCs, and one might speculate whether a few of M31 
ECs are higher luminosity analogues of some of the un- 
usual ultra-faint systems that have recently been discovered 
aroun d the Milky Way, such as Segue I (|Belokurov et al.] 
l2007l ). Some of the ultra-faints which most resemble ECs in 
terms of their luminosities and sizes are labelled in Figure 
1101 There is still much debate about the true nature of the 
lowest luminosity ultra-faint systems, with opinions split be- 
tween their (once?) being ge nuine dwarf galaxies as opposed 
to simple star clusters (e.g. [ Siegel. Shetrone. fc Irwinll2008l : 
iNiederste-Ostholt et al.ll2009l ). In-depth imaging and spec- 
troscopic studies of the M31 ECs as well as the ultra-faint 
galaxy population will provide further insight into this ques- 
tion. 



the MW does not possess compact GCs at large galactocen- 
tric radii. There is also a suggestion of bimodality in the 
size distribution of M31 GCs at large radius with few GCs 
having Rj, in the range from 8 to 15 pc. 

The ECs are of p articular interest since, as first noted 
by iHuxor et al.l l|2005h . they lie between classical GCs and 
dwarf spheroidal galaxies in a plot of My vs R^ (Figure 
nop . The true nature of ECs - whether simply star clus- 
ters or d ark-matter dominat ed systems - remains unclear at 
present. ICoUins et al.l (|2009l ) present a measurement of the 
internal velocity dispersion in HEC12 (EC4) which was used 
to derive a mass-to-light ratio of M/L — Q.Tii]^^ Mq/Lq. 



5 DISCUSSION 

We have presented a detailed study of the properties of the 
M31 GC system, with a particular focus on the halo region, 
using the largest sample compiled to date of confirmed outer 
halo GCs. Wherever possible, we have compared the prop- 
erties of the M31 halo GCs to those of their counterparts 
in the MW, often finding marked differences. In particular, 
M31 has far more GCs than the MW (at least by a factor of 
3) and hosts a population of luminous compact GCs at large 
radius (Rproj < 30 kpc) that, aside from NGC 2419, is com- 
pletely absent in the Milky Way. M31's halo GC system is 
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Figure 10. Plot of against R;, of the ECs (filled 

circles) - with the data from Tanvir et al. (submit- 
ted) used where available, shown with a range of low- 
mass steller sys tems in the Local Group, including: M31 
GCs - circl e s llBarmbv et al.l |2007|'| . MW GCs - crosses 
llHarrisI \l99^ . iKoposov et alj |2007|) , GCs in LMC, SMC, For- 
nax &: Sag, dwarf - crossed circles llvan den Bergh fc MackevI 
12004 iHarrij Il996h . the ECs foun d bey o nd the Local 
Grou p - filled stars l lMouhcine et al.l l20ld iDa Costa et al.l 
20091), UCDs - triangles llMieske. Hilker, fc Infantj |2002| : 
De Propris et a l. 2005; Drink water et al.ll2003h. M31 dwarf galax- 



2004 



2ooe 



les - squares llMc Connachie et al. 2008*; 'Zucker_et_aLl 

20071: iHarbeck et al.. 2005: I rwin et al.l l2008: Martin et al.l 

Ibata et al.ll2007l : iMaiewski et al.ll2007l : lMartin et a l. 2009), clas- 
sical MW dwarf galaxies - aster isks jlrwhr^T^a tzidimitriou 
1 19951 : iMcConnachie fc IrwinI |2006|) . the islolated Local Group 
dwar f s Cetus and Tucana - dotted c ircles l lMcConnachie fc IrwinI 
I2OO6I : ISaviane. Held, fc Piottd [T99I) , and the newly found low- 
lumin osity MW dwarf galaxies - stars jBeloloirov et al.l I2OO8I. 
1201(f). for which t he da ta of iMartin. de Jong, fc Rb3 1I2OO8I ') 
and Ide Jong et al.l i20ldi ) are used. Note that Segue 3 from 
iBelokurov et al. has the structural parameters of a GC, 

and so is plotted as a MW GC. 



also considerably larger in physical extent than the MW's, 
with the most remote member currently-kno wn lying at a 
projec ted (3D) radius of 120 kpc (200 kpcl llMackev et aU 
l2010al ). It is unlikely that the difference in the outer halo 
GC populations is due to intrinsic differences in the shape 
of the GC luminosity functions as we have shown that these 
are in good agreement, although the peak values in the two 
systems are slightly offset. The outer halo of the MW is 



inhabited primarily by Palomar-type GCs which are char- 
acterized by low luminosities and diffuse structures. Such 
clusters may also exist in the outer halo of M31 but, with 
typical Mvo > —5, they are difficult to detect at that dis- 
tance. A detailed assessment of the completeness of our M31 
GC catalogue is required before we can determine whether 
M31 also possesses an excess of Palomar-type GCs at large 
radius and conduct a full analysis of the GCLF. 

M31 hosts a population of luminous extended GCs 
which currently have no known counterparts in the MW. 
With Mv < —6 and Kh > 20pc, such objects should 
have been easily detected with the SDSS out to a radius 
of ^ 300 kpc withi n the Galactic halo (see Figure 10 in 
iKoposov et"al]|2008l ). It may be that similar objects exist 
in the MW halo yet are so few in number that they lie in 
parts of the sky not yet surveyed to sufficient depth. Al- 
ternatively, if ECs represent the tail of the size distribution 
of clusters at a given magnitude, then the absence of lumi- 
nous examples in the MW may simply be due to the overall 
difference in the number of luminous GCs in the two sys- 
tems. It is possible that some of the Palomar-type GCs in 
the MW are the low luminosity analogues of the ECs seen in 
M31. It is of interest to note that many of the Palomar-type 
GCs have been shown to have younger ages than the bulk of 
the halo GCs and evidence supports the notion that many 
of these objects formed in dwarf galajcies that su bsequently 
merged with the MW (e.g. lMackev fc van den Bergh 2005 ). 
On the other hand, although EC4 does not appear to pos- 
sess a dark matter halo, it cannot yet be ruled out that some 
ECs in M31 may share more of a kinship with dwarf galaxies 
- perhaps representing the bright tail of the population of 
ultra- faint dwarfs that has recently been uncovered in the 
MW (|Belokurov et al.ll2007l ). Some of the ECs are very faint 
and extended (e.g. HEC13, see Table 1), approaching Will- 
man 1 in its properties. 

The M31 GC population has often been considered to be 
somewhat redder than that of the MW but with the redden- 
ings we have adopted here the two systems have very similar 
colours. The halo GC populations (Rproj ^15 — 30 kpc) in 
the two systems are especially noteworthy as they have al- 
most identical (V-I)o colours which remain constant with 
radius. Using a subsample of M31 GCs with known [Fe/H] 
from CMD-fitting, the mean metallicity of the outer halo 
GCs is [Fe/H]~ —1.9 ± 0.2 dex with no discernible gra- 
dient. This is in reasonable agreement with the mean for 
the outer MW GC population, [Fe/H]~ -1.7 ± 0.2 dex 
beyond 15 kpc, but considerably lower than the mean of 
the M31 stellar halo over the same radial range, [Fe/H] 
-0.8 to —1.4 dex (Kalirai et al. 2006: Cha pman et al.l 
Richardson et al.,,2009. ). This argues 



l2006l : lKoch et al.ll2008l : 



against the outer M31 field halo and halo GC system forming 
in situ at the same epoc h. Fur t herm ore, following the argu- 
ment made by Scarlc fc ZinnI l| 19781 ) in the context of the 
MW, the lack of a metallicity gradient within the M31 halo 
GC population argues against the system having formed as 
part of a pressure-supported slow collapse. Instead, it is con- 
sistent with the picture wherein the outer GC system formed 
in a number of smaller subsystems which later merged to 
form the halo. 

The large radial baseline spanned by our sample of GCs 
has enabled construction of the GC number density profile 
to hitherto unprobed distances. An unexpected result is the 
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flattening of the proflle beyond a radius of 30 kpc (see Fig- 
ures [7] and [8]). Intriguingly, a similar flattening has been 
observed in the surface brightness profile of t he underly- 
ing st ellar halo. This was first pointed out by llrwin et alj 
l|2005h who used RGB star count data in a large swath 
centered on the souther n minor axi s of M 31 and later con- 
firmed and extended bv llbata et al.l l|2007h who constructed 
an azimuthally-averaged stellar surface brightness profile for 
the entire south-east quadrant of the galaxy, reaching dis- 
tances of ~ 150 kpc. Figure [Til compares the M31 GC num- 
ber density profile wi th the metal - poor minor axis profile of 
the stellar halo from llbata et al.) (|2007l ). Both profiles flat- 
ten at roughly the same radius, 25 — 30 kpc, however the 
GC number density proflle appears flatter bey ond this break 
than that of the field stars, llbata et al.l (|2007l ) showed that, 
beyond a projected radius of 30 kpc, the stellar halo proflle 
can be flt with a power-law of form Ey oc ^-i-9i±o i2 ^jiQg 
we have found here that the GCs behave as oc Ji-0-S7±o.52 ^ 
It is unclear at present whether this difference should be 
viewed as signiflcant. The slope of the stellar surface density 
profile is highly dependent on accurate subtraction of con- 
taminating foreground stars, while the GC number density 
profile suffers from small number statistics at large radius. 

The fact that both GC and fleld star popula- 
tions reveal profile fiattenings at around the same ra- 
dius is consistent with the idea that accretion has played 
an i mportan t role in the formatio n of the outer halo. 
lAbad i. Nav arro, fc Steinmet j l|200d ) have used numerical 
simulations of galaxy formation within a ACDM cosmologi- 
cal framework to investigate the structure of galaxies formed 
through a combination of in situ star formation and ac- 
cretion. They show that beyond the luminous edge of the 
galaxy, defined as that radius at which the accreted stel- 
lar component starts to dominate over that of the in-situ 
component, the slope of the radial surface brightness profile 
changes (see their Figure 4). Within this radius, which corre- 
sponds to 20 kpc in their simulations, the surface brightness 
profile is well fit by a de Vaucouleurs bulge plus an expo- 
nential disk. At larger radii, the outer halo profile flattens 
significantly and can be flt with a power-law which varies 
from Ey cx R'^-^ at 30 kpc to oc fl"^" at 100 kpc. This be- 
haviour is qualitatively consistent with both the outer fleld 
star and GC number density proflles in M31, however the 
observed radial fall-offs are shallo wer than the simulations 
predict (see also llbata et al.l (|2007l )V It is intriguing that the 
GC number density proflle in M31 has an outer slope that is 
reasonably close to that expected for the dark matter halo 
- oc R~^'^ - calculat ed over the range 25-100 kpc, usin g the 
parameters given in iKlvpin. Zhao, fc Somerv"iMl2002l ). 

The lack of a metallicity gradient in the halo GC popu- 
lation and the shape of the GC areal number density proflle 
both provide support for accretion playing a signiflcant role 
in building up the M31 GC system. The sheer number of 
outer halo GCs and the existence of particularly extended 
clusters could also be signatures of this mode of formation. 
In this scenario, the overall differences in the halo GC pop- 
ulations of the MW and M31 could be the result of the two 
galaxies having experienced a different number of accretion 
events, or accretions of a different type. For example, the 
MW may have accreted mostly low mass satellites which 
carry few, if any, associated GCs while M31 may have under- 



gone at least one more substantial merger fe.g. lFardal et all 
2008). 

While ample evidence exists for satellite accre- 
tion events contributing fleld stars to stellar ha- 
los 



lllbata. Gflmore. fc Irwin 

iFerguson et al.l |2002| : iMcConnachie et al 



1994 



Ibata et al.l I2OOII : 
20091 ). direct 



evidence for GC systems being built up in this manner 
has been less forthcoming. One notable example is the 
Sagittarius dwarf that is currently being accreted onto the 
Milky Way and which has been shown to be contributing 
at least one massive compact GC (M54) as well as several 
Palomar-typc clusters (Bcllazzini, Forraro, fc Ibata 20031; 
lOa Costa fc Armandroflj 1 19951 : iForbesl^BridgesI l20iol ). 
More recently, direc t evidence for GC accr etion in M31 has 
been presented by iMackev et al.l (l2010bl ). These authors 
examine the spatial correlation between the positions of 
a sample of halo GCs (many of which are included in the 
present sample) and underlying tidal debris streams. They 
use a Monte Carlo approach to show that the probability 
of the observed degree of alignment being due to chance 
is low, below 1%, and conclude that the observed spatial 
coincidence reflects a genuine physical association. They 
further argue that the accretion of cluster-bearing satellite 
galaxies could plausibly account for > 80% of the M31 
GC population beyond 30 kpc. The properties of the M31 
halo GC population presented in this paper are wholly 
consistent with this idea. 

Finally, it is interesting to compare the halo GC sys- 
tem of M31, the most populous and extended GC system of 
a disk galaxy in the local Universe, with that of the well- 
studied giant elliptical M87. Wide-fleld ground-based studies 
of M87's GC populat ion ha ve recently be en carried out by 
iTamura et al] (|2006al lbl) and lHarrisI (|2009| ) enabUng the sys- 
tem to be traced out to distances of > 100 kpc. Since GCs 
at the distance of M87 (~ 16 Mpc) are unresolved from the 
ground, these studies need to rely on statistical subtraction 
to detect and characterize the GC population which c an lead 
to som e un certainties, particularly at large radii. Both lHarrisI 
(2009") and lTamura et al.l (|2006bl ) flnd that the red GCs in 
M87 are concentrated in the central regions of the galaxy 
(<^ 50 kpc) while the blue population can be traced out to 
at least 100 kpc and shows no evidence for a colour gradient 
in these parts. This is very similar to our finding of a large 
colour spread in the inner regions of M31 with a very uni- 
form blue population dominating in the outer halo. The M87 
studies also derive radial number density profiles for the GCs 
and show a single proflle form (either R^^* or power law R~" 
where n smoothly increases with radius) flts the data well 
out to > 150 kpc. However the proflle shape at large radius 
is highly dependent on the assumed level of contaminating 
sources (e.g. foreground stars, compact background galax- 
ies) and it is not yet possible to place rigorous constraints 
on the form of M87's GC distribution in these parts (al- 
though the GC system of M87 has the beneflt of not suffer- 
ing from problems associated with s mall number st a tistics) . 
Indeed, inspection of Figure 5 from iTamura et al.l l|2006bl ) 
hints at a possible flattening in the number density proflle 
at 90 kpc, s imilar to what we see at ~ 30 kpc in M31. 
ITamura et akl (12006^ ) also note that the blue GCs in M87 
are more extended than the stellar light, but r ecent work by 
iJanowiecki et all l|2010D : IWilliams et al.l (|2007l ) shows M87's 
stellar halo is far larger than previously thought. Using very 
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Figure 11. Plot of the M31 GC radial profile (diamonds) and 
the arbitrarily-scaled minor-a xis metal-poor (— 3.0 < [Fe/H] < 
—0.7) profile (asterisks) from llbata et al.l ([SoO^)- Errors for the 
stellar profile are less than, or comparable to, the symbol size. 
The GC data are presented in the same bins as in Fig. [7] The 
flattening toward the central regions of M31 seen in that figure 
was more apparent than here as it employed a logarithmic axis 
for galactocentric radius. 



deep wide-field imagery, they trace the M87 surface bright- 
ness profile to ~ 180 kpc at which point the V-band surface 
brightness falls below their detection limit of 29 magnitudes 
per square arcsec. This is entirely consistent with the ex- 
tent of the blue GC population. In summary, although M87 
and M31 differ vastly in terms of their inner morphologies 
and galaxy environments, their halo GC populations show 
some striking similarities. Perhaps the properties of halo GC 
systems are determined by processes which are largely de- 
coupled from those which shape the main components of 
galaxies as we see them today. 



6 SUMMARY 

We have investigated the global properties of the M31 GC 
system using an updated sample which includes newly- 
discovered GCs reported by us in Paper I as well as other 
revisions to the RBC. We also derive structural parameters 
for 13 ECs. We find that many of these are less lumin o us an d 
less extended than those presented in iHuxor et all (|2005h . 
bridging a gap between them and the "Palomar-type" GCs 
four id in the MW. 

iMackev et al.l (|2010bl) recently showed that many clus- 
ters in the outer regions of M31 are physically associated 
with tidal streams, and the results presented in this paper 
are entirely consistent with this scenario. Specifically, we 
find no evidence for a significant radial colour/metallicity 
gradient at large galactocentric radii, as expected from 
accretion. We also find evidence for a fiattening in the 
GC number density radial profile in M31 occurring at 
a projected radius of ~ 30 kpc, coinciding with a sim- 
ilar feature in the underly i ng st ellar halo component. 
lAbadi. Navarro, fc Steinmet j l|2006l ) have shown that such 
a flattening occurs naturally in galaxies that grow through 



a combination of in situ star formation and accretion, with 
the point of transition indicating the radius beyond which 
the bulk of the matter has been accreted. 

Wherever possible, we have compared the properties of 
the M31 halo GC system to that of the MW, often find- 
ing marked differences. Although the overall form of the 
luminosity functions is similar in both systems down to 
Mvo ~ —5, M31 possesses a significant population of lumi- 
nous and compact GCs at large galactocentric radii which, 
aside from NGC2419, have no counterpart in the MW. M31 
also has a number of extended GCs, many of which are far 
larger than those in the MW (FigurefTO)). On the other hand, 
halo GCs in M31 and the MW have similarity blue mean 
colours beyond Rproj > 15 — 30 kpc, with little dispersion, 
indicating that old metal-poor populations dominate in both 
cases. We suggest that the differences between the two GC 
systems could be, at least partly, explained by the differing 
accretion histories that M31 and MW have experienced. 

Finally, our work illustrates the importance of extend- 
ing study of GC systems to large galactocentric radii. The 
accretion of cluster-bearing satellites is likely the dominant 
process in building up the halo GC populations of galaxies, 
while in situ formation may contribute the most to the pop- 
ulations at smaller radii. As a result, it is the halo GC pop- 
ulations that have the most to tell us about the hierarchical 
assembly of galaxies. Large area surveys with telescopes such 
as Pan-STARRS and LSST should contribute significantly 
to building samples of GC candidates in the remote out- 
skirts of galaxies within and beyond the Local Group, while 
spectroscopic campaigns will become increasingly necessary 
to weed out contaminants in these sparsely-populated parts. 
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